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ABSTRACT 

The radiative efficiency of super-Eddington accreting black holes (BHs) is explored for 
magnetically-arrested disks (MADs), where magnetic flux builds-up to saturation near the 
BH. Our three-dimensional general relativistic radiation magnetohydrodynamic (GRRMHD) 
simulation of a spinning BH (spin a/M = 0.8) accreting at ~ 50 times Eddington shows a to¬ 
tal efficiency ~ 50% when time-averaged and total efficiency > 100% in moments. Magnetic 
compression by the magnetic flux near the rotating BH leads to a thin disk, whose radiation 
escapes via advection by a magnetized wind and via transport through a low-density channel 
created by a Blandford-Znajek (BZ) jet. The BZ efficiency is sub-optimal due to inertial load¬ 
ing of field lines by optically thick radiation, leading to BZ efficiency ~ 40% on the horizon 
and BZ efficiency ~ 5% by r ~ 400 r g (gravitational radii) via absorption by the wind. Im¬ 
portantly, radiation escapes at r ~ 400/',., with efficiency 77 ~ 15% (luminosity L ~ 50Thdd), 
similar to 77 « 12% for a Novikov-Thorne thin disk and beyond 77 < 1% seen in prior GR¬ 
RMHD simulations or slim disk theory. Our simulations show how BH spin, magnetic field, 
and jet mass-loading affect the radiative and jet efficiencies of super-Eddington accretion. 
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1 INTRODUCTION 

Black hole (BH) accretion drives a broad range of phenomena, in¬ 
cluding quasars, active galactic nuclei, BH X-ray binaries, some 
gamma-ray bursts, and tidal disruption events (TDEs). The disk’s 
gravitational potential energy and BH spin energy are converted 
into radiation, winds, and relativistic jets via magnetic stresses gen¬ 
erated by the magneto-rotational instability (MRI) (JBalbusjfcHaw- 
ley jl 998| , magnetic field threading the disk (Blandford & Payne 
1982), or magnetic field threading the BH (Blandford & Znajek 
1977). Magnetic torques are maximized when magnetic flux ad¬ 
verts inward and piles-up until magnetic stresses balance incoming 
gas forces - the so-called magnetically-arrested disk (MAD) state 
(Narayan et al.|2003|. Simulations of non-radiative MADs show up 
to ~ 300% total efficiency (Tchekhovskoy et al.|2011||McKinney| 
|et al.|20 i~2) , corresponding to 3 times the energy out as going into 
the BH. A fundamental question is: What form does this energy 
take at large distances and what fraction is radiative? 

The radiative efficiency of BH systems is expected to be de¬ 
pendent upon the accretion rate M, which controls the density, 
disk thickness, and dynamical importance of radiation (from gas- 
pressure dominated at low M to radiation-pressure dominated at 
high M). To scale Me 1 or the radiative luminosity L, for a black 
hole mass M, speed of light c, gravitational constant G (giving 
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gravitational radius r g = GM/c 2 ), and Thomson electron scatter¬ 
ing opacity K ts , one can use the Eddington luminosity 
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One can also choose to normalize M by M^aa = (l/77NT)TEdd/c 2 , 
where t/nt is the nominal accretion efficiency for the Novikov- 
Thorne thin disk solution ( [Novikov & Thorne[1973| ) (commonly, 
a fixed t/nt = 0.1 is used, but we include the spin dependence). 

For luminosities L > 0.3L E dd> the accretion flow is expected to 
become geometrically thick and optically thick, and in this regime 
the photons can remain trapped within the flow as in the “slim 
disk” model (which includes no magnetic field) (Abramowicz et al.| 


1988). The super-Eddington accretion regime where M > M^i or 


L <: Gtdd may help explain ultra-luminous X-ray sources as highly 
super-Eddington stellar-mass BHs (Watarai et al. 2001, Miller & 
|Colbert|[2004] l. Also, some black hole X-ray binaries can have 
L > L Edd (e.g„ SS433, |Margon et al.| 19791 |Takeuchi et al.|2010[ 
GRS1915+105, [Fender & Belloni 2004k while tida l disruption 


events seem to require M » A?Edd (Bloom et al. 


2011 


In this paper, we seek to test whether super-Eddington MADs 
around spinning BHs are radiatively efficient. We use the gen¬ 
eral relativistic radiation magnetohydrodynamic (GRRMHD) code 
HARMRAD, which uses the M1 closure for radiation (McKinney) 
|eTaLl|20l4) . The slim disk model and recent radiative GRMHD 
simulations of non-MAD (i.e. lower magnetic flux than MAD on 
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the BH and in the disk) super-Eddington flows are radiatively inef¬ 
ficient with L < L em even for M ~ \0QM Eii (S^dowski et al.|2014 


[McKinney et al.j[2014j >, far below the thin disk Novikov-Thorne 
(NT) efficiency (e.g. about 6 % efficient for BH spin of a/M = 0 
to 12% efficient for a/M = 0.8). However, MADs (not yet studied 
with radiation) can maximize the efficiency and may lead to a ra¬ 
diatively efficient state. This is also plausible for rapidly spinning 
BHs because MADs are then compressed into a thin disk l }McKin-| 
ney et al. 20121 , which can lose radiation more rapidly. 

We discuss the physical and numerical setup in ^2] Results 
and discussions are presented in ^3] We summarize in SR] 


2 FULLY 3D GR RADIATIVE MHD MAD MODEL 

This model with black hole mass of M = 10M o and dimensionless 
spin a/M = 0.8 has nominal thin disk efficiency ~ 12.2%, so 
that Me<m ~ 1-2 x 10 19 g/s (see Eq.[IJ. While we choose a specific 
BH mass, the flow is strongly electron scattering dominated and so 
the results might apply roughly equally to all BH masses. 


2.1 Opacities 

We assume solar abundances (mass fractions of Hydrogen, Helium, 
and “metals”, respectively, X = 0.7, Y = 0.28, Z = 0.02) giving 
electron fraction Y e = (1 + 30/2 and mean molecular weight jj. x 
0.62, entering gas entropy, pressure, and temperature T g [Kelvin], 
We use a frequency (v) mean of the opacity a v to get an 
absorption-mean absorption opacity (units of [cm 2 /g]) of k = 
(f dva v J v )/(f dvj v ), where the absorbed radiation J v is assumed 
to be a Planck distribution at a radiation temperature of T r = 
(E/a m , i ) 1/4 [Kelvin], where E = (/'w v R /JV , iC is fluid 4-velocity, R 
is radiation stress-energy tensor, and a m d is the radiation constant. 
The electron scattering opacity is 

K CS * * 0.2(1 + X) Kkn , (2) 


where the Klein-Nishina correction for thermal electrons is K kn 
(1 + (T g /(4.5 x lO 8 )) 086 ) -1 ( jBuchler & Yueh|l976 1 . 

The absorption-mean energy absorption opacity is 

Kabs ~ (- 1 -) • (3) 

\ Km + KH- ^Chianti + Kbf + Kff J 


Buchler & Yueh 1976 


which bridges between different temperature regimes. The molec¬ 
ular opacity is K m ~ 0.1 Z, and the H~ opacity is kh- ~ 1.1 x 
10~ 25 Z°- 5 p 0 5 rJ' 7 . Rest-mass density, p, is in cgs units. The Chi¬ 
anti opacity accounts for bound-free at slightly lower temperatures 
and is given by K C hianti ~ 4.0 x 1 0 34 p(Z/Z mhT )Y l .T g ]1 T; 30 (this ac¬ 
counts for the assumed Z = Z so]ar = 0.02 for figure 34.1 in |Draine] 
(2011 ', most applicable for baryon densities of n b ~ 1cm -3 ). The 
bound-free opacity is 


K b{ ~ 3 x 10 25 Z(1 +X + 0.75TjpTg 05 77 30 In (l + 1.6 (T r /T g )), (4) 

( jRybicki & Lightman|1986| >, where the 1 + 3C+0.75 Y term is roughly 

accurate near solar abundances, and the ln() term comes from fitting 
the absorption-mean integral. The free-free opacity is 


Kff 


* 4 x 10 22 (1 +20(1 - Z)pT g 0 ' 5 T/ 3 ' 0 In (1 + 1.6 (T r /T g )) 

x (l +4.4x 10~' 0 r g ), (5) 


for thermal electrons and no pairs (see Eq.5.25 in |Rybicki & Light^] 
|man|1 986 and Shu| 1 99l). Thermal energy Comptonization is in¬ 
cluded as in Spdowski et al.| ([2015|>. The mean emission opacity 


Kemit is the same as /c abs but letting T r —> T g , such that Kirchoff’s 
law gives an energy density emission rate of A = cpK e ^ t a IiA T g . The 
total opacity is K tot = K es +/r a b s . The low-temperature opacities avoid 
unphysical opacity divergences during the simulation. 

2.2 Initial Conditions 

The initial disk is Keplerian with a rest-mass density that is Gaus¬ 
sian in angle with a height-to-radius ratio of H/R » 0.2 and radi¬ 
ally follows a power-law of p oc r -0 - 6 . The solution near and in¬ 
side the inner-most stable circular orbit (ISCO) is not an equilib¬ 
rium, so near the ISCO the solution is tapered to a smaller density 
(p —> p(r/15) 7 , within r = 15r s ) and a smaller thickness (H/R —> 
0.2(r/10)°' 5 , within r = 10 r g - based upon a low-resolution simu¬ 
lation). The total internal energy density u tot is estimated from ver¬ 
tical equilibrium of H/R =3 c s /vk for sound speed c s ss -y/T tot P tot /p 
with T tot 4/3 and Keplerian speed v K « (r/r g )/((r/r g ) 3/2 + a/M). 
The total ideal pressure P tot = (T tot - l)M tot is randomly perturbed by 
10% to seed the MRI. The disk gas has T gas = 5/3. The disk has an 
atmosphere with p = 10 -5 (r/r s ) -1 1 and gas internal energy density 
£gas = 10 ~ 6 (r/r g y 5/1 . The disk’s radiation energy density and flux 
are set by local thermal equilibrium (LTE) and flux-limited diffu¬ 
sion ([McKinney et al.|2014J with a negligible radiation atmosphere. 

We do not use polish donuts ( |Abramowicz et al.|1978j > or equi¬ 
librium tori as initial conditions. The outer parts of tori have an ex¬ 
tended column of gas at high latitudes that contributes significantly 
to spurious luminosity (see section 6.8 in |McKinney et al.|2014) l, 
trapping of the disk wind and radiation, and artificial (instead of the 
self-consistent wind) collimation of the jet. 

The initial magnetic field is large-scale and poloidal. For r < 
300r s , the coordinate basis (/(-component of the vector potential is 

A 0 = MAX(r'T0 40 - 0.02,0)(sin d) l+h , (6) 

with v = 1 and h = 4. For r > r 0 = 300r g , the field transitions 
to monpolar using A$ = MAX(r/jl0 40 - O.O2,O)(sin0) 1+ ' ,(, '° /r) . The 
field is normalized with ~ 1 MRI wavelength per half-height H giv¬ 
ing a ratio of average gas+radiation pressure to average magnetic 
pressure of /3 ~ 24 for r < 100r s . 

2.3 Numerical Grid and Density Floors 

The numerical grid mapping equations and boundary conditions 
used here are identical to that given in |McKinney et al.] l [2012| >. The 
radial grid of N r = 256 cells spans from R in » 0.688r H (horizon 
radius jh) to R out = 10 5 r s with cell size increasing exponentially 
till r ~ 500r g and then even faster. Radial boundaries use absorb¬ 
ing conditions. The 0-grid of N e = 128 cells spans from 0 to n 
with mapping given in [McKinney et al~[ ( |2012| ) but with n ]a = 0.7 
to follow the jet, and other coefficients are slightly tuned so that 
the grid aspect ratio at r ~ 30 r g is 1 : 2 : 3. Approaching the 
horizon (where the disk thins due to magnetic compression), the 
grid is tuned to follow the compressed disk, such that on the hori¬ 
zon there are 20 points across a half-height of the actual final disk 
with thickness of order H/R ~ 0.1. The Poynting-dominated po¬ 
lar jet contains, respectively, about 80,60,110 0 grid cells near the 
horizon, r ~ 20 r g , and r > 500r g . This gives sufficient resolution 
of the Poynting-dominated jet. The 0-grid of N$ = 64 cells spans 
uniformly from 0 to 2n with periodic boundary conditions. 

As in our other papers i jMcKinney et al.|20T2||2013||2014} , we 
test the so-called convergence quality factors for the MRI in the 6 
and (/> directions (<2 o,mri and S^.mri) and turbulence (Q n i m?cor ) mea¬ 
suring, respectively, the number of grid cells per MRI wavelength 
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in the 9 and <p directions and the number of grid cells per correlation 
length in the radial, 9 and (/> directions. At late times, our simulation 
has 20 ,mri ~ 340, <2^mri ~ 45, and rest-mass and magnetic energy 
densities have 2„/ m , cor ~ 25,20,6 at r ~ 8 r g , indicating good r, 9 
resolution and marginal (p resolution. We also measure the thick¬ 
ness of the disk per unit MRI wavelength, Sj, where Sj < 0.5 is 
where the MRI is suppressed. Initially ~ 0.6 at all radii, while 
the time-averaged flow has Sj ~ 0 . 1 out to r ~ 60r g . 

The rest-mass and internal energy densities are driven to zero 
near the BH within the jet and near the axis, so we use numerical 
ceilings of b 2 /p = 300, (> 2 /e gas = 10 9 , and e gas /p = 10 10 . The value 
of b 1 jp is at the code’s robustness limit for the chosen resolution. 
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2.4 Diagnostics 


The disk’s geometric half-angular thickness ( H ) per radius (R) is 

8 8 i/W 

with n = 2, H 0 /R = 0, and surface differential dA G( p. One computes 
9q like H/R, but let n = 1 and [9q, Hq/R ) = rr/2. 

The mass accretion rate and energy efficiency are, respec¬ 
tively, 


M 

1 



f(E + pu r + R\)dA g ,p 

W\h 


( 8 ) 

(9) 


where T is the gas stress-energy tensor and [M] H is the time- 
averaged M on the horizon. The specific angular momentum ac¬ 
creted is ] = (f(T£ + R r ^)dA e ^)/[M] H . Both 77 and j are composed 
of free particle (PAKE=kinetic-(-gravitational), thermal (EN), elec¬ 
tromagnetic (EM), and radiation (RAD) components. The jet is 
defined as PAKE+EN+EM, in locations where magnetic energy 
density exceeds rest-mass energy density. The wind is defined as 
PAKE+EN, located outside the jet. The wind can also contain un¬ 
tapped EM and RAD components. The dimensionless magnetic 
flux is 


f dA e ,p0.5\B'\ 
T * 0.7-- ~ 


( 10 ) 


for field B in Heaviside-Lorentz units (McKinney et al. |2014j l. 

To obtain the optical depth, at each instant in time we compute 

T* f P« m dl. (11) 

For the radial direction, dl = -f y dr,f y ~ u'(l-(v/c)cos9),(v/c) ~ 
1 - l/(w ') 2 (as valid at large radii), 9 = 0 , and the integral is from 
r 0 = 3000 (being some radius beyond which only transient material 
would contribute to the optical depth, but a radius the disk wind 
has reached) to r to obtain r r (r). For the angular direction, dl = 
f y rd9, 9 = jt/ 2, and the integral is from each polar axis toward the 
equator to obtain r g (9). The flow’s “true radiative photosphere” is 
defined as when t, = 1 , a conservative upper limit to the radius of 
the photosphere for an observer, because radiation can escape by 
tracking with relativistic low-density parts of the jet. 

The radiative luminosity is computed at each instant as 

L = ~f dA e* R t’ (12) 

which is usually measured at r = 400r g , where we only include 
those angles where the gas has 7 >(r) < 1 . 


Figure 1. The initial (t = 0) state consists of a weakly magnetized radially- 
extended vertically-Gaussian disk with H/R ~ 0.2 around a spinning 
(a/M = 0.8) BH. Rest-mass density is shown as color with legend, while 
green lines show magnetic field lines. Density values are made dimension¬ 
less using MEdd. length r g , and time r g /c. The initial disk is threaded by 
weak (but ordered) magnetic flux, capable of accumulating onto the BH 
and leading to the MAD state. 


3 RESULTS AND DISCUSSION 

Fig.[T]and Fig.[2]show the initial and final state of the accretion flow. 
The initial disk is threaded by a weak large-scale poloidal magnetic 
field around a spinning black hole with a/M = 0.8. Rotation am¬ 
plifies the magnetic field via the MRI leading to accretion of mass, 
energy, angular momentum, and magnetic flux. The magnetic flux 
accumulates and eventually forms a quasi-steady super-Eddington 
(M ~ 400LEdd/c 2 ~ 50AfEdd) state that is a MAD (where the MRI 
is suppressed) out to r ~ 60r g after the model is evolved for a time 
31,200 r g /c. The magnetically-compressed radiatively efficient thin 
disk near the BH is exposed by the jet channel, and the magnetized 
wind carries a significant fraction of radiation away from the disk. 

Fig. [3] shows the mass accretion rate, efficiencies, disk thick¬ 
ness, and magnetic flux as time-averaged from t = 30, 000r g /c to 
31,200 r g /c. The disk is in inflow and thermal equilibrium out to 
r ~ 20r g with constant fluxes of mass, energy, and specific angu¬ 
lar momentum vs. radius. The total (kinetic + gravitational + elec¬ 
tromagnetic + radiative) efficiency r] x 50%. As in non-radiative 
MAD models fTchekhovskoy et al.|201 l[[McKinney et al.|2012| >, 
the efficiency is beyond the Novikov-Thome value of t/nt ~ 12%. 

The magnetic field threading the BH is strong, with T « 10 at 
late time and T« 8 on average with B z oc r ~ 5/4 through the disk. 
The magnetic field threading the BH and disk leads to a magnetized 
wind that carries a significant amount of radiation away from the 
disk and helps to avoid the classical photon trapping effect of slim 
disks ( [Abramowicz et al.| 1988fr . The magnetized wind at r ~ 400r g 
has a matter efficiency of about 15%, but the optically thick wind 
also contains about 15% trapped radiation energy and 10% electro¬ 
magnetic energy that could be tapped for acceleration or heating. 
By r ~ 400r g , the wind carries most of its radiation within a half¬ 
opening angle of 30° around the polar axes. 

The Blandford-Znajek (BZ) effect leads to an electromagnetic, 
BZ, efficiency 77 ~ 35% on the BH. Some of that energy forms a jet 
with 77 ~ 10% at r = 507" g , but some of that energy is absorbed by 
the wind, leading to a jet having only 77 ~ 5% by r = 400r g . 

Radiative-loading of magnetic field lines threading the BH 
leads to a lower-than-optimal BZ efficiency (relatedly, see |Beskm] 
|et al.|2004||Takahashi & Qhsuga|201 5| once the MAD state builds- 
up. The inertial loading of magnetic field lines is due to a total 
energy density p tos = p + e gas + E once r > 1 caused by p keep¬ 
ing up with the MAD’s higher b 2 when the numerical ceiling of 
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Figure 2. The evolved (t = 31, 200r g /c) super-Eddington MAD state, with 
the figure split at x = 0 as two panels (x < 0 on left and x > 0 on right). Left 
panel shows radiation-frame radiation energy density (color, with legend), 
radiation velocity lines (black, fixed line width), and optical depth of unity 
away from each polar axis, Tg = 1 (cyan lines). Right panel shows fluid- 
frame rest-mass density (color, same legend), magnetic field lines (black, 
thicker lines for more magnetized gas), where magnetic energy is equal to 
rest-mass energy density (red lines), and same optical depth (cyan lines). 
The MAD state reaches a quasi-steady state out to r ~ 20r g . Radiation 
advects inward within the equatorial disk, outward through the low-density 
jet channel, and outward along with the optically thick wind. 


b 2 /p = 300 is enforced. The optically thick radiation slows the jet 
magnetic field line rotation rate down by order unity compared to 
the optimal BZ value, because fr/Ptot ~ 1 in the funnel (even if 
b 2 Ip » 1 there) as determined by the rough condition that the jet 
and disk have similar b 2 , yet Y » 1 implies that b 2 ~ p in the disk, 
and radiative energy density E ~ p (from a thermal estimate of 
H/R ~ cJvk ~ 0.3). A restarted simulation with an exponentially 
reduced (as b 2 Ip approached its ceiling) opacity has a BZ efficiency 
o f T) ~ 100% on the BH, as consistent with Y » 8 (jTchekhovskoy| 


et al. 


2011 


. A restarted simulation with a lower ceiling b^Jp = 100 
led to a 20% BZ efficiency. So, mass-loading and opacity physics 
in the funnel are important to BZ-driven jets and winds for super- 
Eddington flows where radiation crosses magnetic field lines. 

Fig. 0 shows the jet, wind, and “true” photosphere at large 
radii. The jet-wind boundary oscillates in angle to order unity, caus¬ 
ing the wind to absorb some of the jet and radiation energy. We have 
resolved the true photosphere at r ~ 400r s , and the radiation escap¬ 
ing to an observer has a high radiative efficiency of 77 » 15% (com¬ 
parable with the NT value of 77 ~ 12%), corresponding to luminos¬ 
ity L ~ 50.LEdd- Most of this radiation is within a half-opening an¬ 
gle of 15° around the polar axes. The luminosity in this GRRMHD 
MAD simulation is much higher than predicted by the slim disk so¬ 
lution at 77 ~ 2% or seen in other viscous non-MAD or non-MAD 


MHD simulations at 77 < 1% 1 

Ohsuga et al. 2005J Ohsuga & Mi- 

neshige 2011; S^dowski et al. 

2014J McKinney et al. 2014). Far 

beyond the photosphere, radiation might be absorbed (Sadowski 


|& Nar ayan 2015) or some might leak out of the wind containing 
77 » 15% in trapped radiation, but our grid is too unresolved to 
accurately track energy conversion for r > 600r s . The radiation 
ultimately released is high in our simulation because 1 ) the disk 
near the BH is forced to become thin; 2) a magnetized wind pulls 
radiation-filled material off the disk; 3) the jet drills a channel and 



10 100 
r tr,] 

Figure 3. The time-averaged angle-integrated mass flux, efficiency, disk 
thickness, and magnetic flux. From top to bottom, panels are: Total mass 
accretion rate (M) per unit M E dd out to where there is inflow and thermal 
equilibrium (r ~ 20 r g , beyond which the line is truncated), energy effi¬ 
ciency r\ (total: solid black line, electromagnetic: blue dashed, matter with¬ 
out rest-mass (i.e. kinetic+gravitational+thermal): dark red dotted, radiation 
(negative within 3 r g ): red long-dashed, radiative luminosity L: magenta dot- 
dashed), disk thickness H/R (evolved state: black solid, initial state: blue 
dashed), and magnetic flux (total on BH and threading disk: black solid, 
only threading disk: dark red dotted). The disk is magnetically compressed 
near the BH, and the thin disk generates a significant radiative flux that ef¬ 
fectively releases starting at r ~ 30()/\,. The initial disk with H/R ~ 0.2 has 
puffed up to H/R - 0.3 from r ~ 20-70 r g after evolving for several thermal 
times. By r ~ 400r s , the luminosity reaches L ~ 50Z,Edd- 


pushes back the opaque disk wind; and 4) the jet becomes conical 
at large radii allowing free-streaming of radiation. These lead to an 
effective diffusion timescale shorter than in the slim disk model. 

We also performed otherwise identical simulations that are 
MAD with a/M = 0 as well as non-MAD simulations using a 
toroidal field in the initial disk (with y3 ~ 20 for all radii at the equa¬ 
tor). The non-spinning MAD never produced a jet and the radiation 
remains trapped by the broad disk wind. Similarly, the toroidal field 
models with a/M = 0 and a/M = 0.8 show no low-density funnel 
region. Similar to |Sadowski & Narayan| l |2015| , our a/M = 0 MAD 
model with M ~ 20M Edd has a wind with NT-level efficiency of 
order 5% and negligible radiation flux at large radii by r ~ 400 r g 
within the photosphere, because all radiation energy flux converted 
to kinetic energy flux. Our toroidal field models with a/M = 0 
and a/M = 0.8 have M ~ Meaa with a wind efficiency of 2% 
and 4%, respectively, and a radiative efficiency of about 1%. So 
the formation of a low-density jet channel by a spinning BH, with 
enough magnetic flux to launch a jet and strong wind, helps super- 
Eddington accretion become radiatively efficient. 

Convergence testing was performed by restarting each 
previously-mentioned model at a resolution twice lower in each di¬ 
mension. The restarted simulation starts at half-way through the 
higher resolution simulation, and then the two resolutions are com¬ 
pared for the latter half of the simulations. We find that the total, 
BH, jet and radiative efficiencies, and Y agree to within 30% or 
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Figure 4. The evolved (t = 31,200 r g lc) super-Eddington MAD state at 
large radii showing the jet, wind, and photosphere. Shows fluid-frame rest- 
mass density (color, with legend), where magnetic energy equal to rest-mass 
energy density (red lines), radiative photosphere at r r = 1 (black line), and 
the wide-angle wind component’s magnetic field line, within which (toward 
the polar axis) the outflow has achieved more than a single flow time along 
the line extending ~ 1000r g (blue line). Quantities are 0-averaged. The jet 
extends to r ~ 30,000r g with Lorentz factor y ~ 5 by r ~ 10 3 r g . The wind 
extends to r ~ 3000r^, with most of the wind’s trapped radiation having an 
opening angle of 30° at r ~ 400/g. The disk at R » 100r g has not evolved 
much. Radiation within r < 400 r g can more freely stream within the jet’s 
low-density channel. By t ~ 30,000 r g /c, at r ~ 400r g within 15° around 
the polar axis, the outflow become conical and optically thin. 

smaller error. Our convergence quality factors for the MRI and tur¬ 
bulence also suggest the simulations are converged. 

A non-relativistic MHD simulation with a/M = 0 by |Jiang| 
|et al,| < |2014] ) measured a NT-level radiative efficiency of about 
5% that they attributed to MRI-driven magnetic buoyancy verti¬ 
cal transport of radiation. Their simulation had a small vertical ex¬ 
tent ±60r g , so r computed from Eq. ED does not include the ex¬ 
tended wind. Over extended distances, a significant portion of ra¬ 
diation energy flux can convert into wind energy flux (|Sadowski &| 
Narayan 2015 h In our a/M = 0.8 MAD simulation, we attribute the 
high radiative efficiency to magnetic compression of the disk into a 
thin MAD, magnetized wind advection of radiation away from the 
disk, and the formation of a low-density jet channel through which 
radiation can more freely stream. Our simulation boundary is at 
r ~ 10 5 r s , with accurate energy conversion out to r ~ 600/ s . Our 
photosphere at r ~ 400r s accounts for wind material that reached 
r > 300(bg. If we ignore the extended wind, for the a/M = 0.8 
MAD model we would miscalculate the radiative efficiency to be 
77 ~ 30% - an overestimate by a factor of two. For the a/M = 0 
MAD model, measuring at r ~ 60 r g gives 77 ~ 5% - an overesti¬ 
mate from rj ~ 0 % (caused the wind absorbing radiation energy). 


4 SUMMARY 

We have performed fully 3D simulations of super-Eddington ac¬ 
cretion, including a simulation with M ~ 400L Edd /c 2 ~ 50M Edd 
onto a rotating black hole with a/M = 0.8. Sufficient magnetic flux 
was distributed throughout the disk that the BH and disk reached 
MAD levels, where magnetic forces pushing out balance gas forces 
pushing in. The MAD state enabled the super-Eddington accretion 
flow to reach its maximum efficiency, with the total efficiency mea¬ 
sured to be about 50% (higher for higher BH spins). Importantly, 


the system has a high radiative efficiency of about 15% (luminos¬ 
ity L ~ 50L Edd ) beyond the resolved photosphere at r ~ 400r s . 
This occurs because the magnetized wind carries radiation away 
from the disk. Also, the magnetized jet creates a low-density chan¬ 
nel for radiation to more freely stream and pushes away the more 
opaque wind. These effects increase the radiative flux escaping the 
disk and diminish conversion of radiation energy flux into kinetic 
energy flux of the wind. This mechanism allows high radiative ef¬ 
ficiencies from super-Eddington accretion systems. 

As applied to jetted TDEs (Tchekhovskoy et aL]|2014} , our 
super-Eddington MAD simulations show how the jet and radiative 
efficiencies depend not only upon BH spin and magnetic flux but 
also upon jet mass-loading physics that leads to sub-optimal BZ 
efficiencies via dragging of BH field lines by optically thick radia¬ 
tion. Also, while the wind absorbs some jet energy, the jet channel 
exposes the (otherwise obscured) inner hot X-ray emitting disk. 
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